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NONCONFORMING FINITE ELEMENT METHODS ON 
QUADRILATERAL MESHES 
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^ . Abstract. It is well-known that it is comparatively difficult to design noncon- 

forming finite elements on quadrilateral meshes by using Gauss-Legendre points 
on each edge of triangulations. One reason lies in that these degrees of freedom as- 
^O . sociated to these Gauss-Legendre points are not all linearly independent for usual 

expected polynomial spaces, which explains why only several lower order noncon- 
forming quadrilateral finite elements can be found in literature. The present pa- 
■^^ , per proposes two families of nonconforming finite elements of any odd order and 

[^ ' one family of nonconforming finite elements of any even order on quadrilateral 

(-H ■ meshes. Degrees of freedom are given for these elements, which are proved to be 

well-defined for their corresponding shape function spaces in a unifying way. These 
elements generalize three lower order nonconforming finite elements on quadrilat- 
erals to any order. In addition, these nonconforming finite element spaces are 
shown to be full spaces which is somehow not discussed for nonconforming finite 
►^ elements in literature before. 

in 

o 
en 

^ I Because of their flexibility and stability when compared with conforming finite ele- 

ment methods, nonconforming finite element methods have become very important 
and effective discretization methods for numerically solving, among others, high 

k> ' order elliptic problems, Stokes-like problems and Reissner-Mindlin plate bending 

H ■ problems. 

" " " Quadrilateral meshes are very important in scientific and engineering computing. 

Indeed, many popular softwares for computations of the fluid mechanics in two 
dimensions are defined on quadrilateral meshes. However, most of nonconforming 
finite element methods for second order problems are defined on triangles [H 131 IH |5l 
[7] while there are only a few nonconforming finite element methods on quadrilaterals. 
Compared with nonconforming triangular finite elements, it is more difficult to 
construct nonconforming quadrilateral finite elements. In fact, a sufficient condition 
for convergence of consistency error terms is to require nonconforming functions to 
be continuous at Gauss-Legendre points on interior edges of the triangulation used 
[T8| [T9| |20| I2T] . Hence, for m order nonconforming finite elements, there are 4m 
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Gauss- Legendre points on the boundary of each element. However, these degrees of 
freedom on Am Gauss-Legendre points are not all linearly independent for the space 
of polynomials whose restrictions on four edges are polynomials of degree < m, see 
( 12. ip for details. This means that at least one relation holds, which motivates the 
use of higher order monomials of one variable in shape function spaces. In such 
a spirit, a class of nonconforming quadrilateral elements is proposed in literature, 
which includes, the Han element [8J, the nonconforming rotated Qi element due to 
Rannacher and Turek [16], the Douglas-Santos-Sheen- Ye (DSSY) element [6], the 
enriched nonconforming rotated Qi element due to Lin, Tobiska and Zhou [13]. All 
of these nonconforming quadrilateral finite element methods are of first order and 
are stable for the Stokes problem with pure Dirichlet boundary conditions for the 
velocity. However, since the usual Korn inequality does not hold for them, they 
can not be applied to the Stokes problem with mixed boundary conditions for the 
velocity. Moreover, these nonconforming elements are somehow not extended to 
higher order nonconforming elements in literature so far . 

In [15], a nonconforming linear element is introduced on quadrilateral meshes, 
which is motivated by a key observation that any linear function on a quadrilateral 
can be uniquely determined at any three of the four midpoints of edges, which leads 
to a set of nodal basis functions whose values at four edges satisfies the aforemen- 
tioned relation. A similar element (but different on general quadrilateral meshes) is 
designed in [9] , which is based on an observation that a frame of the linear function 
space on the reference element can be mapped and glued to form a basis of the 
interpolated space of the conforming bilinear element space by the canonical inter- 
polation operator of the nonconforming rotated Qi element [16]. This explains why 
the element therein is named as "constrained quadrilateral nonconforming rotated 
Qi element" by its authors. Very recently, the ideas of [SI [15] are extended to design 
a nonconforming cubic element on quadrilaterals in [H] . We refer interested readers 
to [ini |12] for the nonconforming quadratic element on quadrilaterals, which can be 
regarded as the quadrilateral counterpart of the triangular Fortin-Soulie element [7]. 

The purpose paper is to propose nonconforming finite elements of any order on 
quadrilaterals in a unifying way, which is somehow missed in literature. Compared 
with nonconforming finite elements on triangles [H IH Hj O [7] , there are at least two 
more difficulties on quadrilaterals: (1) what shape function spaces should be used for 
nonconforming finite elements of any order; (2) how to prove unisolvency of these 
degrees of freedom which consist of values of polynomials at aforementioned 4m 
Gauss-Legendre points and other degrees of freedom in the interior of elements. It 
should be stressed that nonconforming finite elements on quadrilaterals in literature 
are defined and analyzed one by one, see[6l|8l|9l|10l[T2l[T3l[Tl[l5l[l6]. 

In order to overcome the first difficulty, we propose to use the following two 
families of shape function spaces 

(1.1) Rm{K) := PmiK) + span{£'"y - xf"}. 



(1.2) ER^{K) := P^{K) + span{x™y - xy™ , x"'^^ - ^+'1, 

for odd integer m > 1. Here and throughout this paper, Pm{M) denotes the space 
of polynomials of degree < m over the domain M; Qm{M) denotes the space of 
polynomials of degree < m in each variable. For even m, we propose to use the 
same shape function spaces as the serendipity elements of [21 H], namely, 

(1.3) Ri{K) := P^K) + span{x"y, xy™ }• 

These three families of elements generalize the Pi nonconforming element of [9l 
[T5] and the nonconforming cubic element of [H], the nonconforming rotated Qi 
element of [IE] , and the nonconforming quadratic element of [IHl US] to any order, 
respectively. 

Degrees of freedom are given for these shape function spaces, which are proved 
to be well defined based on some observation concerning Legendre polynomials on 
four edges. In addition, these finite element spaces are proved to be full spaces. 

The rest of the paper is organized as follows. Sections 2, 3 and 4 present three 
families of shape function spaces and their corresponding degrees of freedom, which 
are proved to be well-defined. Section 5 defines nonconforming finite element spaces 
which are proved to be full spaces, and shows approximations of these spaces. Section 
6 analyzes consistency errors, which is followed by numerical examples in the final 
section. 

2. The first family of Shape function spaces 

Since Gauss-Legendre points on each edge will be used to define degrees of freedom 
and consequently continuity for new elements under consideration, given an integer 
k > 0, let gi, i = —k, ■ ■ ■ , —1, 0, 1, ■ ■ ■ ,k, denote zeros of Legendre polynomials of 
degrees 2fc + 1 on the interval [—1, 1]. By skew symmetry, Qi = —Q-i, i = 1,- ■ ■ ,k, 
and go = 0. 

Before we present shape function spaces on K := [—1,1]^ and corresponding 
degrees of freedom, we investigate a special relation of values at 8k + 4 Gauss- 
Legendre points Gr,i = (1,5'i), Gi,i = (-1,5'i), Gt,i = (giA), Gb,i = (fi'i,-l), i = 
—k, ■ ■ ■ ,k on four edges of K, of polynomials over K whose restrictions on four 
edges are polynomials of degree < 2k + 1. 

Lemma 2.1. Let v be a polynomial over K such that its restrictions on four edges 
are polynomials of degree <2k + l. Then it holds that 

(2.1) ^^iU{l,gi) + v{-l,gi)\ = ^7J{)((7i,l) + f)((7i,-l)j, 

i=-k ^ ^ i=-k ^ ^ 



where 








(2.2) 


7i = 


k _, fc+1 -, 

1 ^i 1 - - ^ 9j _ 

i^j=-ik+l) ^' ~ ^1 i^3=-k ^' ~ 3^ 


_2 ^ 1-g] 
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for i = —k, ■ ■ ■ , —1, 1, ■ ■ ■ ,k, and 

(2.3) 7o = 4n^- 

Remark 2.2. For k=0, the relation (12. 11) is the constraint used in [51 [U]; for k=l, 
the relation (12. ip recovers that of [14j . 

Proof. The main idea is: Since the restrictions of v on four edges are polynomials of 
degree < 2/c + 1, they can be exactly represented by 2A; + 2 Lagrange interpolation 
basis functions of degree < 2A; + 1 in one dimension, which gives two representations 
of values of v at each corner of /T, then the desired result follows. In order to 
accomplish this, first add the point Q-k-i = — 1 to these 2k + 1 Gauss-Legendre 
points, and define a first set of Lagrange interpolation basis functions: 

k 

£«(£)= Yl ^^^^ioTxe[-l,l],t = ~ik + l),-k,--- ,k. 
i^j=-{k+i) ^' ^^ 

Then, add the point gk+i = 1 to these 2A; + 1 Gauss-Legendre points, and define 
another set of Lagrange interpolation basis functions: 

fc+i 

£i(x) = TT i^—^ioTxe[-l,l],i = -k,---,k + l. 
,4JL, 9^ - 9j 

Since span{/:^(^^^)(x), ■ ■ ■ ,£°(x)} = span{£i;.(x), • ■ ■ ,Ci+i{x)} = P2fe+i([-l, 1]), 
and the restriction of v on each edge of if is a polynomial of degree < 2A; + 1, these 
restrictions can be expressed as 



v{±l,y)= Yl ^'iy)n±'^^9^), v{£,±l)= Yl 'C°(£)^(^.,±l), 

i=-(k+l) i=-{k+l) 

^ ' ' k+1 k+l 

v{±l,y) = Y C]{y)vi±l,g,), v{x,±l) = Y 'Cl{x)v{g„±l). 

i=—k i=—k 

Let ai = jC°(1), i = —{k + 1), ■ ■ ■ , fc, and /3j = C]{—1), i = —k, ■ ■ ■ ,k + 1. The skew 
symmetry of Gauss-Legendre points, and definitions of a_(fe+i) and f3k+i, yield 



(2.5) -.«, - n ^ - n 1^ - A 



j=-k ^^ j=-k ^^ 



Since values -0(1, 1) (resp. v{—l, 1), '0(1, —1) and v{—l, —1)) of two representations 
are identical, it follows from (12 .4^ that 

k k k+l k 

2^ aiv{l,gi)= J^ aiv{gi,l), ^l3iv{l,gi)= J^ aiv{gi,-l), 
i=-(fc+i) j=_(fc+i) i=-k i=-{k+i) 

k k+l k+l k+l 

^ aiv{-l, gi) = ^ I3iv{gi, 1), ^ At'(-1, gi) = ^ Pivigi, -!)• 

j=— (fc+l) i=—k i=k i=—k 

By (I2.5p . this gives 

A: fc 

^{ai + l3i){v{l,gi)+v{-l,g,)) = ^ («i + A)({)((7„ 1) + t;(^„ -1)). 

i=~k i=—k 

Then a direct calculation completes the proof. D 

By canceling a common factor, f l2.ip can be simplified slightly that 

k , 

i=-A; ll|JMi=l(5'i 5' J 

where 



^j 1 , / o 



2 if 2 = 0, 

Ijgl ifz^O. 

Given an odd integer m = 2A; + 1 > 0, recall shape function spaces: 

(2.6) Rm{k) := P„(i^) + span{x"^y - xiT} for any (x, y) G ET := [-1, 1]' 

Note that for fc = 0, Rm{K) is the shape function space of p], i.e., 

i?i = span{l,a;,y}. 



For /c = 1, Rm{K) is the shape function space of 

To define degrees of freedom for this space, let G denote the set of Gauss-Legendre 
points on four edges of K, namely, 

(2.7) G := {(1, g,), (-1, g,), {g,, 1), {g,, -l),t = -k,---, k}. 
To define other degrees of freedom in the interior of K, let 

(2.8) / := {(£,, m)J = 1, ■ ■ ■ , (2A; - 1)(A; - 1)} 

be a set of interior points of K, the standard Lagrange points inside the reference 
element K, so that any polynomial q{x,y) G P2k-z{K) can be uniquely defined by 
its values at points in I. 

Theorem 2.3. If v{x,y) G Rm{K) vanishes at all points in G U I , cf. (12. 7p and 
dMD, thenv{x,y) = 0. 
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Proof. The function v{x,y) G Rm{K) can be expressed as 

(2.9) vix,y) = aoix''+'y-xf''^')+ Yl ^m^'^'- 

i+j<2k+l 

As v{l,y) = at 2A; + 1 Gauss- Legendre points, v{l,y) is a multiple of the (2A;-|-l)-st 
Legendre polynomial, '0(1, y) = CoL2k+i{y) for some constant Cq. We will show that 
the constant cq is zero. In fact, by the continuity of v{x,y) at four vertexes of the 
square. 



■i,y) 



V[X 



-1) 



-Co^2fc+l(y) 
-CoL2k+l{x) 



(2.10) 
(2.11) 

This indicates that 

' Co,o + C2,o + C4,o + C6,o H h C2fc,0 

Co,l + C2,l + C4,l + C6,l H h C2fc,l 

Co,2 + C4,2 + C6,2 + C8,2 H h C2fc-2,2 



v{l,y) for y G [-1,1], 
■0(1, x) for X G [—1, 1]. 



< 



Co,2fc-2 + C2,2fc-2 

Co,2fc-l + C2,2fc-1 

Co,2fc 
, Co,2fc+l 
' Co,o + Co,2 + Co,4 + Co,6 + " 

Cl,0 + Cl,2 + Ci,4 + Cifi + ■ 

C2,0 + C2,4 + C2,6 + C2,8 + ' 



■ + Co,2fc 

■ + Cl,2fc 

■ + C2,2A:-2 



0, 
0, 
0, 

0, 
0, 
0, 
0. 
0, 
0, 
0, 



C2fc-2,2 + C2k-2,2 
C2fc-1,0 + C2fc-1,2 
C2k,0 
, C2k+1,0 

Since co,2fc+i = 0, by comparing the coefficients of y^fc+i j,^ f|2.10p . oq is a multiple of 
the leading coefficient of the Legendre polynomial L2k+i'- 

(4fc + 2)! 



0, 
0, 
0, 
0. 



for j = 0) 
for j = 1) 
for J = 2) 

for J = 2k- 2) 
for J = 2k- 1) 
for j = 2k) 
for J = 2k + 1) 
for i = 0) 
for i = 1) 
for i = 2) 



for i = 2k-2) 
for i = 2A; - 1) 
for i = 2k) 
for 2 = 2A; + 1) 



Co: 



'22fc+i(2A;+l)! 
Again, comparing the coefficients of x'^''^^ in (12. lip 

(4fc + 2)! 



-ao- 



Co 



22fc+i(2fc + l)! 



ao. 



Hence cq = 0. Thus v{x,y) vanishes on the whole boundary dK and can be expressed 
as 

v{x,y) = b{x,y)q{x,y) with q{x,y) G Pzfc-al^), 
where the bubble function b{x,y) = (1 — x^)(l — y^). Finally, since v{x,y) vanishes 
at points in /, q{x, y) = 0. This completes the proof. D 

The above theorem implies that any v{x, y) G Rm{K) can be uniquely determined 
by its values at points in G U /. Though the number of points in G U / is {2k + 
3)(A; + l) + 2, which is 1 greater than the dimension {2k + 3){k + l) + l of Rm{K), the 
relation (12.11) implies that the number of linearly independent functionals defined 
for the space Rm{K) is equal to the dimension of the shape function space. This 
motivates the following degrees of freedom for the shape function space Rm{K): 

• values at points in G which satisfy the relation (12. ip ; 

• values at points in /. 

Remark 2.4. We can take the following shape function spaces 

(2.12) Rm{K) := RmiK) + span{x"'y}, and Rm{k) := P^{k) + span{xy'^}. 

3. The second family of shape function spaces 

As we see in the previous section, functionals defined by values at points in G U / 
are not all linearly independent for shape function spaces Rm{K) defined in the 
previous section, there is a relation (12. ip for all functions in Rm{K). To make these 
functionals linearly independent for some shape function spaces, we propose to use 
higher order monomials of one variable, say a;2'^+2 and y^^^^, which motivates to 
enrich Rm{K) by span{x'^^~^'^ — y'^^'^'^}. This leads to the following shape function 
spaces: 

(3.1) ER^k) := PUk) + spanjx^'^+iy - xf'^'^x'''-' - f'^'}. 

Remark 3.1. For k = 0, the space ERm{K) is the shape function space of the 
nonconforming rotated Qi element from [16] .' 

ERi{K) = span{l, x, y, x^ — jf}. 

For the space ERm{K), the degrees of freedom are: 

• values at points in G, defined in (12. 7p : 

• values at points in /, defined in (12. 8p . 

Theorem 3.2. Ifv{x,y) G ERm{K) vanishes at points in G U I , then v{x,y) = 0. 
Proof. By the definition of ERm{K), v{x,y) can be expressed as 

(3.2) v{x,y)=Mx,y) + c,x''^' + c^f'^' + c, {x^'+'y - xf'^') + c,{x''^' - f'^'), 

where vi{x,y) G Pm{K)\span{x'^'^^^ , y^^^^} , and q, i = 1, ■ ■ ■ ,4, are four interpo- 
lation parameters. This implies that the restrictions of v on four edges of K are 
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polynomials of degree < 2k + 2. Since v{x,y) vanishes at m = 2k + l Gauss- Legendre 
points on four edges of K, these restrictions on four edges can be written as 

v{x,l) = L2k+i{x){aix + bi), v{x,-l) = L2k+iix){a2X + b2), 

vi^,y) = L2k+iiy)ia3y + b3), t'(-l,y) = 1^2fc+i(y)(a4y + &4), 

where Oj, 6j, i = 1, ■ ■ ■ ,4, are some constants which will be shown to be zero next. 
To this end, let c = {4k + 2)!/(2^'^"'"^(2A; + 1)!) be the coefficient of the monomial 
^2A:+i q£ L2k+i{x). A comparison of coefficients from f l3.2p and (13. 3p for monomials 

^2fc+2^ ^2fc+2^ X^fc+l^ g^j^^ ^2fc+l Igg^fjg ^Q 

(3.4) ai = 02 = —as = —04 = C4/C, 

and 

Ci + C3 = cfei, Ci-C3 = c62, 
C2- 03 = 063, C2 + 03 = 064. 



(3.5) 



Note that v{x,l) and t)(l,y) are equal at corner (1,1), v{x,l) and v{—l,y) are 
equal at corner (—1,1), v{x,—l) and -0(1, y) are equal at corner (1,-1), v{x,—l) 
and v{—l,y) are equal at corner (—1,-1). Since L2fc+i(l) = — L2fc+i(— 1), these 
observations give 

«i + 61 = 03 + 63, «! — 61 = 04 + 64, 

(32 + 62 = % - 63, a2 - 62 = 0-4 - O4. 

It follows from ([331), (ESI) and (JSJl) that 

(3.7) ai = a2 = 03 = 04 = 61 = 62 = 63 = 64 = oi = 02 = 03 = 0. 

This implies that 04 = and v{x,y) vanishes on the boundary of K. Hence, 

v{x, y) = b{x, y)q{x, y) with q{x, y) G P2A:-3(^), 

where the bubble function b{x,y) = (1 — x^)(l — y"^). Finally, since v{x,y) vanishes 
at points in /, q{x,y) = 0. This completes the proof. D 

For the space ERm{K), we define another set of degrees of freedom as follows 

• moments of order <2k on each edge of K; 

• values at points in /. 

We note that these two sets of degrees of freedom are not equivalent as {2k + 2) 
polynomials are involved in ERm{K). 

Theorem 3.3. For any v{x,y) G ERm{K), suppose its moments of order < 2k on 
each edge of K and values at points in I vanish. Then v{x,y) = 0. 

Proof. By the definition of ERm{K), v{x, y) can be expressed as 

v{x,y) = v^{x,y) + c^x^'^^' + 02^^^=+^ + c,{x''^'y - xf'^') + c,{x''^' - f'-''), 



where Vi{x,y) G Pm{K)\span{x^^^^ , y^^^^} , and q, i = 1, ■ ■ ■ ,4, are four interpola- 
tion parameters. Consider the restriction on edge 62 of -0, denoted by w, which is a 
polynomial of degree < 2k + 2. The function w can be decomposed as 



where 



w = Wi +W2, 



fc+i 

x2i+l 



wi = '^d2ix'^\ and W2 = '^d2i+ix'^ 

i=0 i=0 

where di, i = 0,1,- ■■ , 2{k + 1), are interpolation constants for w. From degrees of 
freedom it follows 

1 /.I 



wix ^dx = 0, i = 0, ■ ■ ■ ,k, and / W2X dx = 0,i = 0,- ■ ■ ,k — 1. 
1 J-i 

Since t&i (resp. ^2) is an even (odd) function on [—1, 1], it holds that 
1 pi 

2«+l^^ _ n „• _ n _ _ _ U c.^A / „?,_^2i 



■u)ia; (ia; = 0, i = 0, ■ ■ ■ , /c, and / tf)2a^ ^dx = 0,i = 0,- ■ ■ ,k. 
-1 i-i 

Note that the degree of polynomial wi (resp. ^2) is not more than 2k + 2 (resp. 
2k + 1). Therefore both wi and W2 are Legendre polynomials (up to multiplication 
constants), namely, 

wi = a2L2k+2{x), and W2 = 62-^2fe+i(£), 

for two constants 02 and 62. Similar arguments apply to the restrictions on the other 
three edges of K, which leads to 

V\ei = 0'lL2k+2{y) + biL2k+l{y), 
V\e2 = a2L2k+2{x) + &2-^2fc+l(£), 

v\e3 = a^L2k+2{y) + hL2k+i{y), 

V\ei= 04-^2^+2 (£) + WL2k+l{.x). 

Since the coefficient before monomial x^^"*"^ is opposite to that before monomial 
^2A;+2 £qj, ^^ -(.j-^jg gives 

(3.8) ai = 03 = —a2 = —04. 

Let c 7^ be the coefficient of monomial x^^+i of L2k+i{x). A comparison of coef- 
ficients for monomials x^^^^, and y'^^^^ in -Olgj, i = 1, ■ ■ ■ ,4, and those in -0, leads 
to 

Ci + C3 = C64, Ci - C3 = C62, 

Note that ■i)(a;, 1) and '0(1, y) are equal at corner (1,1), {)(£, 1) and {)(— l,y) are 
equal at corner (—1,1), v{x,—l) and "0(1, y) are equal at corner (1,-1), {"(x, — 1) 
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and v{—l,y) are equal at corner (—1,-1). Since e = L2fc+i(l) = — L2fc+i( — 1) and 
d = L2fc+2(1) = -^2fc+2(-l), this leads to 

dai + ebi = da^ — eb^, da^ + eb^ = da^ + eb^, 
(3.10) 

da^ — eb^ = da2ecb2, dai — ebi = da2 — eb2. 

It follows from ([3SD, (O and ( IXTOj) that 

ai = a2 = a^ = a^ = bi = b2 = b-i = b^ = ci = C2 = Cs = 0. 

This implies that C4 = and v{x,y) vanishes on the boundary of K. Hence, 

v{x,y) = b{x,y)q{x,y) with q{x,y) E Pafc-sl^), 

where the bubble function b{x,y) = (1 — x^)(l — y'^). Finally, since v{x,y) vanishes 
at points of J, q{x,y) = 0. This completes the proof. 

D 

4. The third family of shape function spaces 

For nonconforming elements of even order, there is always a discrete bubble func- 
tion which vanishes at Am Gauss-Legendre points. Thus, in additional to one extra 
term for Rm of odd m, we need another extra term enriching P^ polynomials so 
that the finite element function has exactly Am degrees of freedom on the element 
boundary. For m = 2k, we define the third family of nonconforming elements by 

(4.1) Rt,{K) := P^{K) + span{x"y,xr}- 

For the continuity requirement, we need use the m = 2k Gauss-Legendre points, 

9-k, ■ ■ ■ ,9-1,91, ■ ■ ■,9k- 

Lemma 4.1. Let v G R^ he a polynomial over K such that its restrictions on four 
edges are polynomials of degree < 2k. Then it holds that 

.^ 2) y^ t)(l, 9i) - v{-'^, 9i) - v{9^, 1) + v{9i, -1) ^ Q 

o^-k 9^i^ - 9l)]\l\^j=i{9^i - 9''j) 

Proof. We add the point 9-k-i = — 1 to the 2k Gauss-Legendre points, and define a 
first set of Lagrange interpolation basis functions: 

k . _ 

A°(^)= n ^^' ^=-{k+i),-k,---,-i,i,---,k. 

0,i^j=-{k+l) ^' * 

At the other end, adding a point 9k+i = 1, we define another set of Lagrange 
interpolation basis functions: 

k+l 

C\[x)= J] ^^, , = -A;,...,-l,l,...,A; + l. 
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Since 

P2fc([-1, 1]) = span{/:°(x), t = -k-l,--- ,-l,l,--- ,k} 
= span{£|(x), i = —k, ■ ■ ■ , —1, 1, ■ ■ ■ , fc + 1}, 

and the restriction of v on each edge oi K is a polynomial of degree < 2k, we have 
the four equations in (12 ■4p . Again, let 



n„Wt^ If « = -*-!, 



a- = £°(1) = ^ ^^"("^) 110,J7^j = -fc gi-gj 



k 1-gj 

1-1 rr^ -^-9j 



9,-1 ll0,i^i = -fc g,-g, II « ^ /t + i, 

Note that, because the Gauss-Legendre points are symmetric, 

a-k-i = Pk+i, and «» = -/3j if i = -k, ■ ■■ ,-1,1,. . .,k. 
By the continuity of v at the four corner vertexes, it yields 

k k 

a-A;-i^^(-l,-l) + 5Z «i^(5'i, -1) = /3fc+i^(l, 1) + X^ A^(l,5'i), 

Oj^i=-k Oj^i=-k 

k k 

a_fc_i{)(-l,l) + ^ ai{)(5fi,l) = a_fc_i{)(l,-l) + ^ ai{)(l,5fi), 

f)^i=-k Oj^i=-k 

k k 

l3k+iv{l,l)+ ^ /3i{)(5fi, 1) = a_fe_it)(-l, -1) + ^ a^v{-l,gi), 

0^i=-k 0^i=-k 

k k 

/3fc+i{)(l, -1) + Y^ (5iv{gi,-l) = Pk+iv{-l,l)+ Y^ p.iv{-l,gi). 

0^i=-k Q^i=-k 

Eliminating four corner values of -0, we get 

k 

Y "^(^(I'fi'i) -^(-l>fl'i) -^{9iA) + v{gi,-l)) = 0. 

0^j=-fc 

This is simplified to ( I4.2p . D 

To define degrees of freedom for this space, let the set of even Gauss-Legendre 
points on four edges of K be 

G+ := {(1, gi), (-1, gi), {g„ 1), {g„ -1), ^ z = -k,- ■ ■ ,k}. 
To define other degrees of freedom in the interior of K, let 

I+:={{xe,yi),i=l,--- ,{2k~3){k-l)} 
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be a set of interior points of K, the standard Lagrange points inside the reference 
element K, so that any polynomial q{x, y) G P2k-i{K) can be uniquely defined by 
its values at points in /+. 

Theorem 4.2. Ifv{x,y) G R^{K) vanishes at all points m G+UJ+U {(1, 1)}, then 
v{x,y) =0. 

Proof. The function v(x, y) G R^{K) is a P2fc polynomial when restricted to y = 1. 
On the edge y = 1, v vanishes at 2k Gauss- Legendre points plus a corner point 
{(1,1)}. SoC' = Oony = l. Repeating the argument, as v is continuous at the four 
corners, we find -0 = on the whole boundary. 

v{x,y) = b{x,y)q{x,y) with q{x,y) G P2fc-4(^), 

where b{x, y) = (1 — x^)(l — y^). Here with a careful division, we find the coefficients 
for x"^y and y"^x in v are zero. So g(x, y) G P2k-i{K). Finally, since ■0(x, y) vanishes 
at points of /, q{x,y) = 0. This completes the proof. D 

The above theorem implies that any v{x,y) G R^{K) can be uniquely determined 
by its values at points in G^ U /^ U {(1, 1)}. But the number of points is one more 
than the dimension of R^{K). The relation f l4.2p implies that the number of linearly 
independent functionals defined for the space R'^{K) is equal to the dimension of 
the space. This motivates the following degrees of freedom for the shape function 
space R'^{K): 

• values at points in G'^ which satisfy the relation (14. 2p : 

• value at point (1, 1); 

• values at points in /+. 

Note that the value at point (1, 1) is to determine the coefficient of the discrete 
bubble function in v: 



k{x,y) = \{{x^ + f-l-g^,). 



i=l 

5. Nonconforming finite element spaces 
This section defines nonconforming quadrilateral element spaces. 

5.1. Quadrilateral Mesh. Let Th '■= {Ki,i = 1, ■ ■ ■ , A^e} be a shape regular 
quadrilateral partition of Q with diam(i^j) < h. We assume that the partition 
Th satisfies the bisection condition of [17]: The distance dx between the midpoints 
of two diagonals of each element K is of order 0(/i^). 

For a given element K E Th, its four nodes are denoted by Ai{xi, yi),i = 1,- ■ ■ ,4 
in the counterclockwise order. Let K := [—1,1]^ denote the reference element with 
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^K 




Figure 1. The reference element K and a quadrilateral element K. 

nodes Aj(xj, ^j), i = 1, ■ ■ ■ ,4, shown in Figured] Define the bilinear transformation 
Fk-K ^Khy 



X 



^XiNi{x,y), y = ^yiNi{x,y), {x,y) e K, 

i=l 1=1 

where Ni{x, y),i = 1,2,3,4 are the bilinear basis functions, which can be written as 

Ni{x,y) 



i(l -x)il-y), N^ix, y) = \{l + x)(l - y). 



Ns{x,y) = -{l + x){l + y), N4{x,y) = -{1 - x){l + y). 

5.2. Nonconforming finite element spaces and dimensions. For an odd inte- 
ger m = 2k + 1 > 0, define nonconforming finite element spaces by 



(5.1) 



Rh := {v G L'^{Q),v\k o Fj^^ G Rm{K) for any K G Th, v is continuous 
at m Gauss-Legendre points of each interior edge of 7^}- 



(5.2) 



ER^ := {v G L^{Vi),v\K ° Fj^^ G ERm{K) for any K G Th,v is continuous 
at m Gauss-Legendre points of each interior edge of T/i}. 



(5.3) 



ER^ := {v G L\a),v\K o F^' G ER^{K) for any K e %, [v]qds = 

J e 

for any q G P2k{e) for each interior edge e of 7/i}, 



where [f ] denotes the jump of v across edge e. We note again that ER^ is different 
from ER^ due to the higher order polynomial term (a;^^"*"^ — y'^^^'^). For an even 
integer m = 2k, the nonconforming finite element space is defined by 



(5.4) 



^h '■= {fJ ^ L'^i^),v\K o Fj^^ G R'^{K) for any K E Th,v is continuous 



at m Gauss-Legendre points of each interior edge oi Th}- 
The corresponding homogeneous spaces are defined, respectively. 
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(5.5) Rhfi '■ = {v & Rh,v vanishes at m Gauss-Legendre points 

of each boundary edge e oi Th}- 

(5.6) ER^Q := {v E ERJ^,v vanishes at m Gauss-Legendre points 

of each boundary edge e oi Th}- 

(5.7) ERf;^'^ := {v G ERff, f vqds = for any q G P2fc(e) 

for each boundary edge e of 7^}. 

(5.8) i?^Q := {v G -R/[,f vanishes at m Gauss-Legendre points 

of each boundary edge e of Th}- 

5.3. Approximations of nonconforming finite element spaces. Given K G 
Th, define 

(5.9) Gk := G o Fk, and Ik ■= I o Fk- 

Then, define the canonical interpolation operator Tl^^^ : H'^{Q) -^ ER^ by 

(5.10) {Yi^^v\K){p) = v\k{p) for any p eGk^Ir and K e % 

for any v G H'^{^). 

Define the canonical interpolation operator Her '■ H'^{^) — > ER^ by 

(n^^i;|x)(p) = v\k{,p) for any p G Ji^, 

for any K e T and t; G iJ^(fi). 

To define an interpolation operator for Rh, let 11*^ be the canonical interpolation 
operator of the conforming Qm element space Qm,h '■= {v G H^{ri),v\K o Fj^^ G 
Qm{K),K G Th}- Then, define an interpolation operator 11^ : H^{Q) — )■ Rh by 

(5.12) {IIrv\k){p) = n'^vlKip) for any p e Gk U Ik a.nd K e % 

for any v G II'^{Q). Since the values at points in Gk of Il'^f satisfy the relation 
02. ip . this operator is well-defined. 
For the even order nonconforming finite elements, we define 

(5.13) iUR+v\K)ip) = n'^vlKip) for any p G G+ U /+ U {F^(l, 1)} and K e % 

for any v G H'^iVt). Since the values of Qm polynomial 11^ v at points in G^^ satisfy 
the constraint (14. 2p . the operator 11^;;+ : H'^{fi) — > R'^ is well-defined. 
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An immediate consequence of these interpolation operators is the following ap- 
proximation property 

(5.14) inf \\Vh{u-Vh)\\Q<Ch"'\\u\\m+i, 

Vh&Vh 

provided that u G H"^^^{Q) and the mesh satisfies the bisection condition where 
Vh = Rh, ER^ , ERf^ , R'^. 

5.4. The full space of nonconforming finite elements. For the finite element 
spaces ER^Q and ERff^, local nodal basis functions are uniquely determined, which 
are glued together to form the global nodal basis functions. Thus, the interpolation 
operators 11^^ and H^ji, cf. fIS.lOp and (15. lip , are on-to mappings, by which the 
dimension of the finite element spaces can be counted. 

Let Ny, Ns, and Ne denote the numbers of vertexes, edges and elements of the 
partition Th, respectively. Let Ny, Ny, Ng and Ng denote the numbers of interior 
vertexes, interior edges, boundary vertexes and boundary edges, respectively. The 
dimensions are 

(5.15) dimE/2f = dimERffQ = NE{2k - l){k - 1) + Ni^{2k + 1). 

However, for nonconforming spaces Rhfi and R'^q, the local nodal basis is not 
unique. The canonical basis functions are not linearly independent but subject to the 
constraints (12. ip and (14. 2p . When gluing these functions together to form the global 
space, it is not clear if they can form a basis for the full space. In other words, the 
local constraints (12. ip (on each element) might be linearly dependent globally, that 
is, some constraints may be automatically satisfied when the neighboring elements 
are subject to the constraints. If so, the interpolation operators 11^^ and 11^^+ in (I5.12p 
and (I5.13P are not on-to mappings. Then, the computational finite element spaces 
are not the full spaces defined mathematically in (15. ip and (15. 4p . but subspaces. This 
problem is not well-addressed in previous research. We give a rigorous analysis. 

Theorem 5.1. The interpolation operators are on-to mappings that 

(5.16) nfi(Q„,/,ni/o'(^)) = ^M, 

(5.17) nH+(Q„,,n/7o'(^)) = ^M- 

The dimensions of these spaces are 

(5.18) dimRh,o = NE{{2k - l){k - 1)) + A^^ + N's{2k), 

(5.19) dimi?+o = NE{{2k - 3){k - 1) + 1) + A^^. + N's{2k - 1). 

Proof. We analyze the space Rhfi only as that for R^q is the same. Rhfi is defined 
as the full space. But the range of the interpolation operator is a subspace 

UR{Qm,hnH^{Q))cRh,o- 
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To show that the two spaces are equal, we show that the dimension of the range of 
Ur is no less than that of the full space Rh,o- The former is easily counted by the 
definition of the operator Ur in f l5.12p . 

(5.20) dimRange(nR) = NE{{2k - l){k - 1)) + Nl^ + N's{2k). 
We define a completely discontinuous space 

Rt = {ve L^in) I v\k o F^^ g RmiK) \/K e %}. 

Then 

(5.21) dim< = Ne ■ (dimP2fc+i + 1) = NE{{2k + 3){k + 1) + 1). 
Now, we introduce linear functionals on the space Rf: 

(5.22) /, : Rt^R\ 

t^ ^„v r,. \v{gd for any QieGx^ dQ, 

(5.23) ji{v) = < , _ 

[i^i9i) -v{g^) for rest g^ G Gk H VL\ 

for any K E Th, where we randomly choose a plus side for each edge, and v{g^) and 
v{g^) are values of v on the two sides of an edge where gi is a Gauss-Legendre point, 
defined in (12. 7p . The nonconforming finite element space Rh is the kernel space of 
the space of above linear functionals. 

Assume the domain is simply connected. Let gi be any one boundary Gauss- 
Legendre point. Let the linear functional evaluating at this point be /q. We show 
next that value of fo{v) = v{gi) is completely determined by the rest nodal values 
for all V G Rh- For example, if the domain consists of one element K, then v{gi) = 
7i"^ X]i=^j ~'yj'^i9j)y cf- (I^-IP - We define a set for the rest linear functionals in (I5.23P 

(5.24) ^ = {/, I ^ = 1, 2, ■ ■ ■ , Nsi2k + 1) - 1}. 



Assume 

(5.25) Yl "^ifi = 0- 



dim J^ 



j=l 



That is 



dim J^ 



J2cJ,ivh) = Vv^eRt. 

i=l 

For convenience, we let cq = 0. Let Vh G Rf and Vh = on all elements except the 
element where /o is defined. We have then 

dimj^ 4m— 1 

= J]^ Cifiivh) = ^ Ci{±Vhigi)). 

i=0 i=0 

By the constraint (12. ip . 

Ci = c-fi or - C7i 
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for some uniform constant c. But above {q} contains cq which is 0. Thus all q 
related to the element K are 0. In this fashion, we know that all q in f l5.25p are 
and the set J-" is linearly independent. Then the dimension of its kernel is dim J-". 
By (lOTD and (KM . (I5l8l) holds, 

dim Rhfi = dim R^ — dim T 

= NE{{2k + 3)(A; + 1) + 1) - {Ns{2k + 1) - 1) 
= NE{{2k -l){k- 1)) + WsW + iV^, 

where we used the relations ANe = 2iV| + N^, 2Ne = 2NI, + N^-2 and N^ = N^^. 
(15161) follows (I5:20|) and (iSlSj) . D 

6. The analysis of consistency errors 

To analyze consistency errors, we need some additional interpolation operators. 
On the interval [—1,1], define two L^ projection operators Vm-i '■ -^^([~1;1]) ~^ 
Pm-j([— 1, 1]), for i = 1,2, respectively. 



(6.1) 



/ Vm-iwqds= wgds for any g G P™_j([-1, l]),i = 1,2. 



for any w e L'^{[-1, 1]). Define the interpolation U^-i ■ C°([-l, 1]) -^ P^_i([-1, 1]) 
by 

(6.2) iUm-iv)igi)=vigi),t = -k,--- ,k, for any t; G C°([-l, 1]), 

where Qi are Gauss-Legendre points on the interval [—1, 1]. We recall that ei and 63 
be two edge of K that parallel to the y axis, see Figure [H 

Since the analysis of consistency errors for both Rh.o and ER^q is similar, only 
details for ER^q are presented as follows. 

Lemma 6.1. Let Hm-i be the interpolation operator defined in (16. 2p . Then it holds 
that 



/ u{v -Iljn-lV\e3)dy - / u{v -Ilm-lV\e:,)dy 
Jes Jei 



< C\^?i\ ■ \n)\ 

— ^ \"'\W^{K)\^\H"^{K)-: 



for any u G H"^{k) and v G ER^{k). 

Proof. Let Vm-i be the L^ projection operator defined in (16. ip . By which, we use 
the following decomposition 

/ u[v-Um-iv\ei)dy- j u[v -Tlm-iv\e^)dy=:Ii + l2 

Je-i Jei 
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where 

l2= Vm-lU\es (^leg " 11^- 1 f) | eg ) rf^ 

-/ Vm-lU\ei{v\ei-'^m-lV\ei)dy. 

The first term Ji can be bounded by the Cauchy-Schwarz inequahty, the trace the- 
orem and the usual Bramble-Hilbert lemma, 

\h\ < C!\u\fjm(^K)\''^\H"^{K)- 

To analyze the second term I2, introduce the following decomposition for v 

V =Vi +V2, 

where -Oi G Pm{K) + span{a;^^+^y — xy'^'^^^}, and V2 = Co(a;^^"'"^ — y2fc+2^ with Cq, an 
interpolation constant. Since -Oilei — nm^i-Oile^, i = 1,3, are polynomials of degree 
at most 2k + 1, and vanish at 2A; + 1 Gauss-Legendre points of e,, 

with two constants q, i = 1,3. Since degrees of polynomials Vm~iu\ei, i = 1,3, are 
not more than 2k, this gives 

(6.3) / Vm-lU\e3{l -'n.m-l)vi\esdy - Vm-iu\e, {l - 'nm-l)vi\e^dy = 0. 

Since -02 lea = "^21(11, it holds that (J — nm„i)-02|e3 = (/ — Ilm-i)v2\ei- Hence, 

/ Vm-lU\e3 {I - '^m~l)v2\e3dy - "Pm-lWlei (-?' " n„_i){;2 |ei# 

(6.4) 

= / Vm-l{u\e3-u\ei) {I -'n^m~l)v2\e3dy. 

From facts that (/ — nm_i)'02|e3 is a polynomial of degree < 2k + 2 and vanishes at 
2k + 1 Gauss-Legendre points of 63, it follows 



(6.5) 



/ Vm-i{u\e3-u\ei) {I -'nm-l)v2\e:idy 
/ {I -Vm~2)Vm^l{u\e3-u\ei) {l - Tlm~l)v2\e3dy ■ 
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It is straightforward to see that the right hand side of (16.51) vanishes for all u G 
Pm-i{K), which leads to 

,^ ^. / 'Pm-l{u\e3-u\ei) {I -'nn,-l)v2\e;dy 

(6.6) J-i 

A summary of (16. 3p . (16. 4p . (16. 5p . and (16. 6p proves 

which completes the proof. D 

Let 62 and 64 be two edge of K that parallel to the x axis, see figure [H A similar 
argument of the above lemma can prove the following result. 

Lemma 6.2. Let Hm-i be the interpolation operator defined in (16. 2p . Then it holds 
that 



/ u{v -Um-lV\e4)dx- u{v - Um-lVleijdx 

J €4, J €2 






for any u E H"^{K) and v E ER^{k). 



Theorem 6.3. Suppose that the quadrilateral mesh Th satisfies the bisection con- 
dition. Then it holds that 

ELtiMds 
(6.7) sup '-^^ . < Ch"^\u\m+i, 

for any u E //'"+H^)- 

Proof. Based on Lemmas 16. II and |6.2[ the proof follows from similar procedures used 
in[9l[n]. D 

Remark 6.4. Similar arguments can show similar estimates for consistency error 
for spaces Rhfl, ERffQ and R^q. 

Theorem 6.5. Let u anduh be the exact solution and finite element solution of the 
Poisson equation, respectively, 

( Vm, Vv) = {f,v) \fvE H'-^' (Q) n H^ (Q) , 
{Vuf„Vv) = {f,v) ^vEVh, 
where \4 = Rh,o, ER^^, ER^^^ or R+q. Then 

\u-Uh\Hi<Ch^'''^'-^"'MHr(n), 
where m = 2k — 1 for Vh = Rh,o, ER^q or ERffQ, and m = 2k for Vh = Rho- 
Proof. It is standard, by applying the Strang lemma, cf. [1]. D 
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7. Numerical test 



We compute the Pk nonconforming finite element solutions on uniform rectangular 
grids for the following Poisson equation: 

-Au = f in fi = (0,1)^ 
u = on dQ. 

In computation, the exact solution is, 

(7.1) u{x,y) = 2\x-x'^){y-y'). 

The first level grid is the unit square. Each subsequent grid is a refinement of the 
last one by dividing each square into 4, denoted by {%,}■ We use the following 6 
nonconforming finite elements: 



(7.2) 
(7.3) 
(7.4) 
(7.5) 
(7.6) 
(7.7) 



K 



{3E) 



K 



(4) 



K 



(5) 



V, 



{5E) 






(6) 



(7) 



{veL^ 
{veL' 
{veL' 

{veL^ 



v\k e Ps 



{xy^ — x^y, x^ 



y'}, He ±P P2}, 



v\k ePi® {xy\x^y}, [v]^ ± P3}, 



v\k e P5 

v\k e P5 



{xy'}, He ± P4}, 

v\KePe®{xy'^,X%}, He ^^5}, 

v\KeP7®{xy'^}, He^^e}, 



P4, 



where K is any square element in Th and e is any edge in the grid Th. To be precise. 



K 



(3E) 



K 



(4) 



K 



(5) 



Rh,0} 



K 



(5E) 






(6) 



(7) 



ER 

p+ 

Rh,Qi 



h,Oy 



defined in (15. 6 p with m = 3, 
defined in (15. 8 p with m = 4, 
defined in (15. 5 p with m = 5, 
defined in (15. 6 p with m = 5, 
defined in (15. 8 p with m = 6, 
defined in (15. 5 p with m = 7. 



We list the computation results in Tables [T]-[6l In all cases except the case of poly- 
nomial degree 3, we stop the computation when the machine accuracy is reached, 
i.e., the relative error of computed solutions is about 10~^^. We can see, for all the 
finite element spaces, the computational solutions converge at the optimal order of 
rate. This is proved in paper. 
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Table 1. The error and order of convergence by the 2D C ^-P^ 
element (OD for ( fTT]) . 





U -Uh\ 1,2 


h' 


\u - Uh Hi 


h' 


1 


0.000000000 


0.0 


0.00000000 


0.0 


2 


0.172089821 


0.0 


1.15734862 


0.0 


3 


0.012510804 


3.8 


0.16038663 


2.9 


4 


0.000823397 


3.9 


0.02155950 


2.9 


5 


0.000052434 


4.0 


0.00280349 


2.9 


6 


0.000003300 


4.0 


0.00035752 


3.0 


7 


0.000000207 


4.0 


0.00004514 


3.0 


8 


0.000000013 


4.0 


0.00000567 


3.0 



Table 2. The error and order of convergence by the 2D C ^-P/^ 
element (I73D for (O). 





U - Uh\ L2 


K- 


U-Uy\Hl 


U 


1 


1.520882827023 


0.0 


10.1100205010 


0.0 


2 


0.073186215065 


4.4 


0.9404045783 


3.4 


3 


0.002467158503 


4.9 


0.0655417961 


3.8 


4 


0.000078057209 


5.0 


0.0042062381 


4.0 


5 


0.000002441049 


5.0 


0.0002645249 


4.0 


6 


0.000000076219 


5.0 


0.0000165552 


4.0 


7 


0.000000002381 


5.0 


0.0000010349 


4.0 



Table 3. The error and order of convergence by the 2D C ^-Ps 
element (EID for f lTl]) . 





1 U -Uh\ L2 


K- 


1^ - W/i Hi 


h' 


1 


0.162340154 


0.0 


0.94025897 


0.0 


2 


0.003652811 


5.5 


0.05082935 


4.2 


3 


0.000061390 


5.9 


0.00176166 


4.9 


4 


0.000000983 


6.0 


0.00005729 


4.9 


5 


0.000000016 


6.0 


0.00000182 


5.0 


6 


0.000000000 


6.0 


0.00000006 


5.0 
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Table 4. The error and order of convergence by the 2D C ^-P^ 
element (OD for ( TTT]) . 





U -Uh\ 1,2 


h' 


\u - Uh Hi 


h"- 


1 


0.160588210 


0.0 


0.93096655 


0.0 


2 


0.003712721 


5.4 


0.05204506 


4.2 


3 


0.000062480 


5.9 


0.00180745 


4.8 


4 


0.000000999 


6.0 


0.00005869 


4.9 


5 


0.000000016 


6.0 


0.00000186 


5.0 


6 


0.000000000 


6.0 


0.00000006 


5.0 



Table 5. The error and order of convergence by the 2D C ^-Pq 
element (ES]) for f lTT]) . 





\\U - Uh\ L2 


h"- 


u-Uh\m 


h"- 


1 


0.051014969254 


0.0 


0.4290029114 


0.0 


2 


0.000428314992 


6.9 


0.0080196383 


5.7 


3 


0.000003402349 


7.0 


0.0001296547 


6.0 


4 


0.000000026614 


7.0 


0.0000020505 


6.0 


5 


0.000000000209 


7.0 


0.0000000322 


6.0 


6 


0.000000000021 


3.3 


0.0000000012 


4.8 



Table 6. The error and order of convergence by the 2D C ^-Pr 
element (IIZD for (TTT]) . 
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